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We report on a search for large extra spatial dimensions in the dielectron and diphoton channels using a data 
sample of 1.05 fb^^ of pp collisions at a center-of-mass energy of 1.96 TeV collected by the DO detector at 
the Fermilab Tevatron Collider. The invariant mass spectrum of the data agrees well with the prediction of the 
standard model. We find 95% C.L. lower limits on the effective Planck scale between 2.1 and 1.3 TeV for 2 to 
7 extra dimensions. 
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Within the standard model (SM) the mass of Higgs boson is 
unstable against radiative corrections. The fact that the mass 
is not of the order of the GUT or Planck scales at 10^^ or 
10^*^ GeV but rather O (10^ GeV) is commonly referred to as 
the "hierarchy problem". One way to circumvent the need for 
such fine tuning in the Higgs mass is by extending the dimen- 
tionality of the space, as in the large extra dimension model 
(LED) proposed by Arkani-Hamed, Dimopoulos and Dvali 



(ADD) 0]], which posits that the fields of the standard model 
are pinned to a (3 + l)-dimensional membrane, while grav- 
ity propagates in nd additional compactified spatial dimen- 
sions. Gauss' Law gives the relation between the fundamental 
Planck scale Mo, the observed Planck scale Afpi, and the size 
of the extra dimensions R: [Mpi] ^ « R''" [Md] '^"'^^ ■ If R is 
large compared to the Planck length ~ 1.6 x lO"'^'' cm, AId 
can be as low as (1 TeV), thus avoiding the hierarchy prob- 
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lem and making gravity strong at the TeV scale. Extra spa- 
tial dimensions will manifest themselves by the presence of 
a series of graviton states, known as a "Kaluza-Klein tower", 
(Grk)- At colliders, large extra dimensions can be probed 
by searchingfor the effect of Gkk on fermion or boson pair 
production [211 . 

Extra dimension amplitudes will result in enhancement of 
the cross sections above the SM values, especially at high en- 
ergies. The LED cross section, which consists of SM, inter- 
ference, and direct gravity terms, can be parametrized by a 
single variable rja = F/M^ where Ms is the the effective 
Planck scale, the ultraviolet cutoff of the sum over Kaluza- 
Klein states in virtual graviton exchange. The exact rela- 
tionship between Ms and Mo depends on the exact quantum 
gravity scenario although they are of the same order of mag- 
nitude. The dimensionless parameter F to leading order (LO) 
and the sub-leading dependence is given by 



F 
F 



1, (GRW d) 
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where s is the center of mass energy of the partonic subpro- 
cess. 

In this Letter, we present a search for LED performed in 
events containing an e+e~ or 77 pair with 1.05 fb~^ of pp 
collider data collected with the upgraded DO detector |5i| be- 
tween October 2002 and February 2006. With 127 pb"i of 
data, DO has published limits on Ms ranging from 0.97 to 
1.44 TeV for Ud = 7 - 2 in the combined dielectron and 
diphoton channels DO has also published limits in the 
dimuon channel with 246 pb~^ of data The efficiency 
and resolution for high energy electromagnetic (EM) objects 
at DO are superior to those for muons and so a search for LED 
in combined e+e~ and 77 (di-EM) final states is superior to 
the dimuon channel. DO and CDF have also published lim- 
its on Mo in the monophoton and monophoton plus monojet 
final states, respectively jH. 

Events are collected using triggers requiring the presence 
of at least one EM calorimeter shower with the transverse 
momentum with respect to the beam axis, px, greater than 
15 GeV. From these data we select e+e~ and 77 events using 
criteria that do not distinguish photons from electrons. We 
require events with two EM showers with px > 25 GeV. 
Showers are labelled CC (EC) if they are reconstructed in 
the central calorimeter (end cap calorimeters) with < 1.1 
(1.5 < \ri\ < 2.4), where pseudorapidity 7^ = — ln[tan(0/2)] 
and 9 is the polar angle measured with respect to the proton 
beam direction. To reduce multijet background, we require at 
least one shower to be in the CC, so that selected events are ei- 
ther CC-CC (both showers in the CC) or CC-EC (one shower 
in the CC and the other in the EC). Each EM shower is re- 
quired to be isolated, with less than 7% of the cluster energy 
in an annular cone 0.2 < ATZ < 0.4 about the shower cen- 
troid, where ATZ = y/ (A?^)^ + (A0)^ and (f> is the azimuthal 
angle. We also demand the scalar sum of the px of all tracks 



in the cone 0.05 < ATZ < 0.4 be less than 2 GeV. Finally, we 
demand the EM shower profile be consistent with that of an 
electron or photon using a test and that 97% of the shower 
energy be contained in the EM calorimeter 

The efficiencies for the electron and the photon selection 
criteria are determined from the same data set used for the 
event selection. We estimate separately the efficiencies for 
the requirement on the EM shower shape, the isolation re- 
quirements based on ATZ, and for all calorimeter-based high- 
Pt triggers relevant to this analysis. In order to estimate the 
different efficiencies, we select a sample of di-EM events sat- 
isfying very loose EM identification requirements with invari- 
ant mass within ±40 GeV around Z boson mass. For each 
of these di-EM candidate events we estimate the efficiency 
as a function of tj versus p^ using the tag and probe method 
10]. This efficiency is then applied to Monte Carlo simulation 
samples. 

The irreducible background to the LED signal is from 
SM e+e" and 77 production, while instrumental background 
arises from multijet and 7 H- jet events with jets misidentified 
as EM objects. To model the invariant mass distribution of 
the physics backgrounds, we use the PYTHIA [ 10] event gen- 
erator using the CTEQ6L1 parton distribution functions ll ill , 
followed by a GEANT-based II12II detector simulation and re- 
construction with the same algorithms as applied to data. The 
next-to-leading order (NLO) effect for both e+e^ and 77 is 
taken into account by multiplying the leading order (LO) cross 
section by a mass independent fc-factor of 1.34 11311 . 

We generate the LED signal for 2 < < 7 and 33 differ- 
ent values of Ms using a parton level generator 11411 . Follow- 



ing 14}], we assume Br (Gkk Br{GKK e+e") = 

2. In order to model the effects of detector resolution and 
initial state radiation (ISR), we generate LEDh-SM and SM- 
only events separately to obtain the parton level distributions 
of the di-EM invariant mass versus the cosine of the scattering 
angle in the centre of mass frame of the two EM candidates 
(I cos 6** I) for each value of Ms and nd considered. The ratio 
of the LEDh-SM and SM distributions are obtained for all val- 
ues of Ms and Ud- Standard model events generated with the 
detailed GEANT-based Monte Carlo simulation are weighted 
by this ratio to model the effect of an LED signal. We reweight 
the shape of the SM to simulate the LED signal, keeping the 
overall normalization as in the pure SM case. By normaliz- 
ing to the Z boson production cross section (NNLO), where 
the signal contribution is negligible, we reduce the fractional 
uncertainty on the product of the efficiency and integrated lu- 
minosity. 

To estimate the normalization of the multijet background, 
we fit the di-EM invariant mass distribution of the selected 
data events with a linear combination of the physics and in- 
strumental background distributions. The shape of the invari- 
ant mass distribution for the instrumental background is esti- 
mated from data events with EM energy clusters that fail the 
requirement for the shower profile. This fit is performed 
in the mass range 60 - 140 GeV where we expect no con- 
tribution from LED. We obtain separate fits for CC-CC and 
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TABLE I: Number of events observed and expected from SM processes in different mass windows for CC-CC and CC-EC events. The 
individual contributions to the total SM expectation from multijet, e^e^ and 77 are also shown separately. 
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FIG. 1: The di-EM invariant mass distributions for CC-CC (a) and CC-EC (b) events. The data are shown by points with error bars, the filled 
histograms represent the Drell-Yan, diphoton and multijet backgrounds, and the solid line represents the total background. The broken lines 
show the invariant mass distributions for two different values of Ms for na — 4. The error bars for the total background include both statistical 
and systematic uncertainties. 
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FIG. 2: The distributions of the center-of-mass scattering angle cos9* of the two final state EM candidates in CC-CC (a) and CC-EC (b) 
events. The data are shown by points with error bars, the filled histogram represent the multijet background, and the solid line represents 
the total background. The broken lines show the cos 6* distributions for two different values of Ms for Ud = 4. The error bars for the total 
background include both statistical and systematic uncertainties. 



CC-EC events. From the fits we determine the fraction /mj of 
the multijet contribution to the total background in the mass 
range 60 - 140 GeV to be /mj = 0.22 ± 0.03 in CC-CC 
events and /mj = 0.24 ± 0.07 in CC-EC events. We ex- 
trapolate the total background using the fitted value of /mj to 
determine the expected number of background events with in- 
variant mass above 140 GeV in both the CC-CC and CC-EC 
configurations. Table |T] shows the numbers of events in dif- 
ferent mass ranges for CC-CC and CC-EC where we would 



expect the LED signal to appear The number of events is 
consistent with the number of expected events from the SM 
expectation. FigurefHa) shows the invariant mass distribution 
for CC-CC events and Fig.IBb) for CC-EC events. The distri- 
butions of I cos 6* I are shown in Fig.|2]for CC-CC and CC-EC 
both for data and the background model. We find that the total 
background distribution for the invariant mass and | cos 6* \ is 
consistent with the data within statistical and systematic un- 
certainties. 
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Most of the systematic uncertainties on the background 
model are dependent on the invariant mass. The dominant un- 
certainty arises from the efficiency of the cut on the shower 
profile used to estimate multijet background ( 1 3% of the back- 
ground itself in CC-CC and 30% in CC-EC). The systematic 
on the LED modeling is dominated by uncertainties on the 
choice of parton distribution functions [(1-19)% in CC-CC 
and (1.5-12)% in CC-EC] . All the other signal uncertainties 
are correlated to SM background uncertainties due to the tech- 
nique used to generate our LED signal. Table HH summarizes 
the dominant background and signal uncertainties taken into 
account in calculating the limit on Mg. The NLO k-f actor un- 
certainty refers to the uncertainty due to choice of PDF, renor- 
malization and factorization scale. 

TABLE II: Systematic uncertainties (in %) on the predicted numbers 
of signal and background events considered in calculating the limit 
on Ms. 



TABLE III: Observed and expected lower limits at the 95% C.L. on 
the effective Planck scale, A4s, in TeV. 
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The two-dimensional distribution of the invariant mass and 
I cos 6* I for the observed dielectron and diphoton events is 
compared with the corresponding distributions expected from 
SM physics and instrumental background, and the LED signal 
for Ms ranging from 1 TeV to 3 TeV for a given Ud- The pos- 
terior probability density P{Ms \ Data) given the number of 
observed events in the fc*'' mass bin and cosO* bin, iV^j,'^', 
is then computed using a Gaussian prior for the SM plus mul- 
tijet background. Evidence of large extra dimensions with a 
given Ms will appear as a peak in P{Ms \ Data) distribution. 
In the absence of signal we proceed to estimate the lower limit 
on Ms using the semi-frequentist CLs method lIlSll . which is 
based on computation of a log likelihood ratio. Both the ex- 
pected and observed limits on Ms at the 95% C.L. are calcu- 
lated. Systematic uncertainties in the signal and background 
distributions are taken into account in this calculation, with 
their correlations properly included. The distribution of the 
ratio of the observed (expected) upper limit at the 95% C.L. 
limit to the predicted cross section as a function of Ms is used 
to extract the observed (expected) limit on Ms for = 7 to 
Ud = 2. 

For the Ud independent GRW formalism, we calculate the 
observed(expected) limit on Ms of 1.62(1.66) TeV. We ob- 
tain the observed limits on Ms at the 95% C.L. in the HLZ 
formalism (sub-leading, rid dependent) ranging from 1.29 to 
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FIG. 3: Observed and expected limits on the effective Planck scale, 
Ms, in the di-EM channel along with previously published limits in 
di-EM channel. 

2.09 TeV for = 7 to = 2. Both the observed and ex- 
pected limits on Ms, for different formalisms and for six dif- 
ferent Ud are summarized in Table Ulll The observed and ex- 
pected limits on Mg for a given number of extra dimensions 
are found to be similar. The present limits are a significant 
improvement over the published limit i^. Figure |3] summa- 
rizes the observed and expected limits on Ms along with the 
previously published limits on Mg in the di-EM channel. 

In summary, we have performed a dedicated search for 
large extra spatial dimensions by looking for effects of virtual 
Kaluza-Klein graviton in the dielectron and diphoton chan- 
nels using 1.05 fb^^ of data collected by DO detector. We 
see no evidence of excess over the standard model prediction 
and set limits at 95% C.L. on the effective Planck scale at 
2.09(1.29) TeV for 2(7) extra dimensions. These are presently 
the most restrictive limits on large extra dimensions. 
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